J2' 

a; 



IFIC/12-67 

Sneutrino Dark Matter in Low-scale Seesaw Scenarios 

Valentina De RomerJ*| and Martin Hirschl 
AHEP Group, Instituto de Fisica Corpuscular - 
C.S.I.C./Universitat de Valencia Edificio de Institutos de Paterna, 
Apartado 22085, E-^6071 Valencia, Spain 

Abstract 

We consider supersymmetric models in which sneutrinos are viable dark matter candidates. 
These are either simple extensions of the Minimal Supersymmetric Standard Model with additional 
singlet superfields, such as the inverse or linear seesaw, or a model with an additional U(l) group. 
All of these models can accomodate the observed small neutrino masses and large mixings. We 
investigate the properties of sneutrinos as dark matter candidates in these scenarios. We check for 
phenomenological bounds, such as correct relic abundance, consistency with direct detection cross 
section limits and laboratory constraints, among others lepton flavour violating (LFV) charged 
lepton decays. While inverse and linear seesaw lead to different results for LFV, both models have 
very similar dark matter phenomenology, consistent with all experimental bounds. The extended 
gauge model shows some additional and peculiar features due to the presence of an extra gauge 
boson Z 1 and an additional light Higgs. Specifically, we point out that for sneutrino LSPs there is 
a strong constraint on the mass of the Z' due to the experimental bounds on the direct detecton 
scattering cross section. 
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I. INTRODUCTION 



Although the existence of dark matter (DM) is supported by a variety of astrophysical 
data, its identity is unknown. Many particle physics candidates have been proposed to 
explain the DM l|. In supersymmetric extensions of the standard model (SM) there are 
the lightest neutralino and the scalar neutrino, which could both provide the correct relic 
density for the DM [2|. The neutralino as a DM candidate has been studied in literally 
hundreds of publications, but also sneutrinos as candidates for the cold dark matter have 

nil 

actually quite a long history |3H5[. However, ordinary left sneutrinos, i.e. the sneutrinos of 
the MSSM (Minimal Supersymmetric Extension of the SM), have been ruled out 6] as the 
dominant component of the dark matter in the galaxy a long time ago due to their large 
direct detection cross section 7|. This leaves only "mixed" sneutrinos, i.e. sneutrinos which 
are partly singlets under the SM group, as good DM candidates. Motivated by neutrino 
oscillation data [8j, in this paper we study scalar neutrinos as DM candidates in models 
with a low-scale seesaw mechanism, either MSSM-like models with an inverse [9] or the 
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11] or models based on an U (1) b-l x U (1)r extension of the MSSM group 



. mear seesaw 

3,Q. 

Singlet sneutrinos as DM have been studied in the literature before. Neutrino masses 
require that pure Dirac sneutrino must have tiny Yukawa couplings. Unless the trilinear 
parameters are hu ge, Dirac (right) sneutrinos are therefore never in thermal equilibrium in 
the early universe 14|, ll5[. 1 However, they could still be non-thermal DM produced in the 
decay of the NLSP ("next-to-lightest supersymmetric particle") [18]. Also, trilinear terms 
are usually thought to be proportional to the associated Yukawa couplings, T u oc Y V A ~ 0(1) 
eV. Treating T u as a free parameter of the order of (9(100) GeV, Dirac sneutrinos can be 



made good thermal DM candidates, as has been discussed in 



19 



sneutrinos of this type have been studied in 



21| . Very light mixed 



221 ] . The LHC phenomenology of mixed Dirac 



sneutrino DM was studied in [23j . Alternatively to a large A-term, Dirac sneutrinos could 



also be made thermal DM in models with an extended gauge group 



24 



25]. 



In the classical seesaw picture 26M29| lepton number is broken at a very large energy 



1 Unless Dirac neutrino masses are due to a tiny vev of a non-standard Higgs field [lfj. In this case, Dirac 
sneutrinos could be the DM and even explain the much discussed claim for a tentative 130 GeV 7 line in 
the FERMI data nj 
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scale, possibly close to the unification scale. In such a setup also the right sneutrinos are 
very heavy and decouple; the sneutrinos remaining in the light spectrum are then very 
MSSM-like. One could, of course, simply put the scale of the seesaw low, say around 
the TeV scale. Yukawa couplings of the order of O(l0~ 6 ) could fit neutrino data and the 
right sneutrinos are thermalized. In such an electro-weak scale seesaw right sneutrinos are 
overabundant unless (i) (again) a large trilinear parameter is assumed [30|; (ii) a new £7(1) 
group is introduced [31]; or (iii) sneutrinos have a large coupling to the NMSSM ("next-to- 
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33]. 



minimal Supersymmetric Standard Model") singlet 

However, the situation is different in extended seesaw schemes like the inverse 
linear seesaw 
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9| or the 



11] . Here, additional singlets need to be introduced, but the neutrino 



Dirac Yukawa couplings can take essentially any value and it is the smallness of the inverse 
or linear seesaw terms which "explains" the smallness of the observed neutrino masses. In 
these setups the sneutrinos are highly mixed states. Inverse seesaw sneutrino DM has been 
studied previously in 34], |35[. Our work differs in several aspects from these earlier papers. 
34j calculated all masses at tree-level and did not carry out a detailed fit to neutrino data, 
while we use full 2-loop RGEs for the parameters, one-loop corrected mass matrices and pay 



35] has some overlap 



special attention to constraints from neutrino masses. Also the paper 
with our work, but concentrates more on collider phenomenology of the inverse seesaw with 
sneutrino DM. 

There are also some recent paper studying extended gauge groups. 



studies inverse 



seesaw in an SU(2)r extension of the MSSM. Also two papers based on sneutrinos in 
^b-l(I) x Uy(1) have been published recently. In [37] an inverse seesaw is implemented 



in Ub-l(X) x Uy(l)- In [38] sneutrino DM within the Ub-l(X) x UyiX) group was studied 
assuming a standard seesaw. However, none of the above papers has studied linear seesaw 
variants. Finally, we mention that part of the results discussed in this paper have been 



presented previously at conferences [c 



All our numerical calculations have been done using SPheno 
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4=11 ] . for which the nec- 



essary subroutines were generated using the package SARAH |42H44] . We have written the 
SARAH input files for the inverse and the linear seesaw, while for the U(1)b-l x U(1)r 
model we used the SARAH input files from The calculation of the relic density of the 
LSP is then done with MicrOmegas 45] version 2.4.5 based on the CalcHep 46] output of 
SARAH. To perform the scans we used a Mathematica package (SSP) 



421. 
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The rest of this paper is organized as follows. In the next section we first recall the 
main features of the supersymmetric inverse and linear seesaws, before discussing briefly 
the minimal U(1)b-l x U(1)r extension of the standard model. In section (lllip we discuss 
phenomenological constraints on the parameter space of the different setups. In section fllVp 
we then calculate the relic density and direct detection cross section. We conclude in section 

(Hi- 
ll. SETUP: LOW SCALE SEESAWS AND EXTENDED GAUGE GROUPS 

In this section we briefly discuss the different setups, which we will use in the numerical 
sections of the paper. We first discuss supersymmetric inverse and linear seesaw, before 
recalling the main features of the minimal U(1)b-l x U(1)r extension of the MSSM. The 
latter can be realized with either inverse or linear seesaw, but has some interesting additional 
features which are not covered by either the inverse or linear seesaw extensions of the MSSM. 

A. Inverse and linear seesaw 

In both, the inverse and the linear seesaws the particle content of the MSSM is extended 
by two types of singlet superfields, v c and S. The former is assigned a L = +1, while the 
latter has formally L = —1. The total superpotential can be written as 

W = W MS SM + W v c + + W hSS (1) 

Here, Wmssm is the usual MSSM superpotential 

Wmssm = Y u uqH u -Y d dqH d -Y e e c lH d +/iH u H d . (2) 

Lepton number conserving terms for the new singlet fields i> c ( "right-handed neutrino" ) and 
S can be written as 

W„c = Y v v c lE u +M R P C S . (3) 

The first term generates Dirac neutrino masses, once the H u acquires a vacuum expectation 
value, while the second term is a mass term for the new singlet fields. In the inverse seesaw 
lepton number is violated by the term 

W ISS = ^ S SS, (4) 
4 



while in the linear seesaw case one writes lepton number violation as: 

W hSS = Y SL SlH u . 



(5) 



In both cases, it is usually assumed that the lepton number violating terms are small 
11] . see also the discussion in section (IHIj) . The neutrino mass matrix and the resulting 
constraints on the model parameters are discussed in section flHI Aj) . 

In supersymmetric models with lepton number violation, also the scalar neutrinos must 



have a lepton number violating mass term 



48] . This term, m 2 M , is given by the difference 



between the eigenvalues of the real and imaginary components of the scalar neutrinos. It is 
therefore convenient to separate the sneutrino mass matrix into CP-even and CP-odd blocks 



49|: 



M" 



(6) 



M\ 
Ml 

Mass matrices for the scalar neutrinos are different in the inverse and linear seesaws. At the 



tree-level, in the inverse seesaw the M\ are given by: 3 



/ 



Mlj S s 



A T 

LR 



m\ + D 2 + {m T D m D ) 

A LR m 2 uC + {M R Ml) + {m D m T D ) 



\ 



m T D M R 
±M R n s + B Mr (7) 
\ M T R m D ±h s MI + B t Mr m 2 +fi 2 +M^M R ±B^ SJ 

Here, D 2 = |m|cos2/9 are the MSSM D-terms, m D = -^v u Y u , A LR = T Yv v u — fim D cotg(3, 
Bm r is the soft bilinear term, Ty v is the soft trilinear and m 2 L , m 2 yC and m 2 s are the scalar 
soft masses for the doublet and the singlets respectively. Only fi s and the corresponding 
bilinear soft term B^ s violate lepton number and only these two come with different signs 
in the CP-even and CP-odd mass matrices. 
For the linear seesaw one finds 



Ml LSS 



( 



A lr ± Ml Ml 



m T D M R ± A T LS 



m 2 L + D 2 + {m T D m D ) + (MjM L ) 

A LR ± M R M L m 2 uC + (M R Ml) + {m D m T D ) ±m D M T L + B Mr 

M T R m D ± A LS ±M L m T D + B t Mr m 2 + M L M T L + M T R M R 



2 Separation into CP-even and CP-odd blocks requires CP-conservation, i.e. all parameters in the mass 
matrices below have to be reaL 

3 We correct some misprints in 34|, |69 1 



with all definitions as in eq. (J7J) and Ml = -j=v u Y SL and A LS = T YsL v u — jj,M L cotg(3. In 
these simple setups all other mass matrices are as in the MSSM and, therefore, not discussed 
here. 



B. Minimal SU(3) C x SU(2) L x U(1) B -l x U(1)r extension of the MSSM 

In order to explain why neutrinos are so much lighter than all other matter particles, we 
have considered in the previous section two variants of the seesaw which can, in principle, 
be implemented at virtually any mass scale. Such seesaw schemes are actually most easily 
realized in a particular class of extensions of the MSSM with an extended gauge group 
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52] based on the SO (10) breaking chains 

SO(10) ->■ SU(3) C x SU(2) L x SU(2) R x U(1) B -l -> SU(3) C x SU(2) L x U(1) Y (9) 
SO(10) ->• SU(3) C x SU(2) L x SU(2) R x U{\) B -l (10) 
->• SU(3) C x SU(2) L x U(1) R x C/(1) B _ L ->• ST/(3) C x ST/(2) L x C/(l) y 

A MSSM-like gauge unification is in this case perfectly viable, and compatible with a U(1)rX 
U(1)b-l stage stretching down to TeV. We will follow eq. (ITU]) , since this variant can be 
realized with the minimal number of additional superfields with respect to the MSSM particle 



content. This model [50l . |52| . which we will call the minimal U(1) b ~l x U(1)r extension 
(mBLR, for short) has been studied in two recent papers . We will follow the notation 



of [13] quite closely. 

The particle content of the mBLR model is given in table (JTJ) . In this setup, the presence 
of v c is required for anomaly cancellation. Breaking the SU (2) L xU(l) b ~l x U(1) b to U (1)q 
requires additional Higgs fields. The vev of the fields xr an d Xr break U(1) B -l x U(1)r, 
while the vevs of H u and break SU(2)l and U(l)y- Note that since H u and are 
charged also under U(1)r, in the mBLR new D-terms are generated in the mass matrix for 
the scalars. These additional contributions with respect to the MSSM allow to have a larger 
mass for the lightest MSSM-like CP-even mass eigenstates and makes it possible to have a 
m h o ~ 125 GeV without constraints on the supersymmetric particle spectrum [12|, 



Assuming matter parity 
the terms 



13], apart from the MSSM superpotential the model also has 



W s = Y u ^LH u + Y s u c xrS - I1rXrXr + UsSS. (11) 





Superfield 


SU(3) C x SU(2) L x U(1) R x U(1) B _ L 


Generations 




Q 


(3,2,0,+i) 


3 




d c 


(3. 1, ~T~77 < ~F ) 


3 




u c 




3 




L 


(1.2.0,-i) 


3 




e c 


(i,i,+i,+i) 


3 




v c 


(l,l,-i,+i) 


3 




S 


(1,1,0,0) 


3 




H u 


(1,2, +|,0) 


1 




H d 


(1,2,-|,0) 


1 




XR 


(1,1, +2. -2) 


1 




XR 


(1.1. -5. +1) 


1 



TABLE I: The Matter and Higgs sector field content of the U(1)r x U(1)b-l model. Generation 
indices have been suppressed. The S superfields are included to generate neutrino masses via the 
inverse seesaw mechanism. Under matter parity, the matter fields are odd while the Higgses are 



even. 



The 2nd term generates Mr = -^Y s v XR while the last term generates the inverse seesaw 
discussed above. The model can, in principle, also be written with a linear seesaw included 
501 ] . Note, that the model assigns lepton number necessarily in a different way then discussed 
in the last subsection, since here B — L is gauged. Thus, B — L is broken by the vevs of \r 
and xr- However, neutrino masses are generated in exactly the same way as in the simpler 
inverse seesaw model, discussed in the previous subsection. 

It is useful to reparametrize the vevs in a notation similar to the MSSM, i.e.: 



•r 



2 1 2 

V XH + V xn 
v Xr 



2 1 2 
Vd+ V u 

Vu 



(12) 



tan/3 R = ^, tan/3 = —. 



J XR 



Vd 



The mass of the new Z'-boson is approximately given by 13 1 

1 



rrir 



{9bl + 9r)v r - 



Z ' <9 2 bl + 9r) ' 4' 
Thus, Vr must be larger than approximately v R > 5 TeV, see also next section. 
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(13) 



Mass matrices for all sfermions for this model can be found in 13j. For us the sneutrino 
mass matrix is most important. In the mBLR model it is given by the expression for the 
inverse seesaw, with exception of M R = -^Y s v XR and new D-term contributions: 

d ll = | {(9bl + 9blr- 9bl9rbl)v 2 r cos(2(3 r ) + (g 2 L + g 2 R + g BL gRBL)v 2 cos2/3) 
d rr = | {[9bl + 9r + 9blr + 9rbl ~ ^9bl9rbl - 2g R g BLR )v 2 R cos(2/3 i? ) 

+(Qr + 9rbl ~ 9bl9rbl - gRgBLR,)v 2 cos 2/3 ) (14) 



Here, D\ L replaces D 2 of the simpler models, while D 2 RR are the new D-terms in the (i> c , z> c ) 
part of the mass matrix. Due to the lower limit for the Z' mass, see eq. ( fl3|) . and since the 
new D-terms in eq. ( 1T4"]) can have either sign, the free parameter tan f3 R is constrained to be 
close to tan/3fj ~ 1, otherwise either one of the sneutrinos (or one of the charged sleptons) 
becomes tachyonic. 

III. PHENOMENOLOGICAL CONSTRAINTS 



In this section we discuss phenomenological constraints on the parameter space of the 
different models. Below, we concentrate on neutrino masses and lepton flavour violation. 
Other constraints on the model space come from SUSY searches at colliders, from Z° phy sics 
(LEP) and from the Higgs results of the LHC collaborations ATLAS 53| and CMS [5^. 

"Heavy" singlet neutrinos with mass below the Z° boson are excluded by LEP experi- 
ments 55|, which set limits on \U"A 2 of the order of to 10~ 5 for the neutrino mass range 



from 3 GeV up to 80 GeV. L3 has searched also for heavy iso-singlet neutrinos decaying 
via N — > IW and set limits which range from \U^\ 2 < 2 x 10~ 3 for masses of 80 GeV to 
Wij\ 2 5-, 10 1 for masses of 200 GeV 56]. Most importantly, the invisible width of the Z° 
boson [57| puts an upper limit on the 3x3 sub-block U^, i,j < 3, of the neutrino mixing 
matrix: 1 - £y= M <i | ELi U ik U Jk \ 



< 0.009 at the 3-cx level even in the case that the 



new mostly singlet neutrinos are heavier than the Z° boson 13] . Finally, the Z° width rules 
out pure left sneutrinos lighter than approximately half of the Z° mass, but sneutrinos with 
suppressed coupling to the Z° below roughly 0.02 — 0.1 with respect to the MSSM coupling 
and masses below < 40 GeV are allowed. 

In inverse seesaw models the Higgs can decay to heavy plus light neutrino, if the heavy 



neutrino has a mass below the Higgs mass 



13 



This limits the Yukawa couplings 



to roughly below \Y V \ < 0.02 for Mr < 120 GeV from measured data on the channel h — > 
WW* — > llvv 58(. For larger masses current Higgs searches provide essentially no constraint 
yet. 

For the model with the extended gauge group searches for a new Z' boson at the LHC 
provide important constraints. Both, CMS 59J and ATLAS [60| have searched for, but not 



observed any hints for, Z n s within the context of different models. For the U(1)b-l x U(1)r 



model the limits are of the order of (roughly^ 



SUSY searches at ATLAS 



6l| and CMS 



m z , > (1.7-1.8) TeV. 



62] provide lower limits on squark and gluino 
masses. For example, in mSUGRA/CMSSM models with tan/3 = 10, A Q = and fi > 0, 
squarks and gluinos of equal mass are excluded for masses below 1500 GeV [61]. This limit 
essentially rules out any value of M 1( / 2 below approximately (600 — 700) GeV for m < 1000 
GeV and My 2 below (350 — 400) GeV in the limit of large mo for pure CMSSM. Of course, 



53 



54j , if interpreted 



the observation of a new resonance with a mass around 125 — 126 GeV 
as the lightest Higgs boson, provides important constraints on SUSY parameters as well. 
However, these constraints are different for the different models we study in this paper. We 
will discuss them therefore when we discuss numerical scans in section (II VI). 



Neutrino masses 



1. Inverse seesaw 



In the inverse seesaw the neutrino mass matrix can be written at tree-level as 



( m T D \ 

m D M R 
\ M T R n s ) 



(15) 



The smallness of the observed neutrino masses is then usually explained as the hierarchy 

US < m D < Mr. 



Following the notation of 63j, we can count the number of physical parameters of the 
model as iVph ys = Ny — Nq + Nq>- Here, Ny is the number of parameters in the Yukawa 
matrices (or mass matrices), G is the original symmetry group which is broken into G' by 
the presence of the Yukawas (or mass terms). In table [Til the counting for the inverse seesaw 
is summarized. 
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Parameters 


Moduli 


Phases 


N Y 


Y e ,Y v ,M R Dirac type 


3 x n 2 


3 x n 2 


N Y 


Us Majorana type 


n(n+l) 
2 


n(n+l) 
2 


G 


U{n) L U{n) v c U{n) e c ® U{n) s 


4 x n(n " 1) 


4 x n(n 2 +1) 


G' 


no LF conservation 






Nphys 




21 
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TABLE II: Parameter counting for MSSM with an inverse seesaw for three generations. 



After absorbing all unphysical parameters by field rotations, we find a total of 30 real 
parameters, 21 moduli (12 masses and 9 mixing angles) plus 9 phases. It is common practice 
to choose a basis in which the charged lepton mass matrix (Yukawa: Y e ) is diagonal, which 
fixes 3 parameters. The remaining parameters could be fixed by going to a basis where 
Mr is real and diagonal. In this case Y v and [is are completely general, arbitrary matrices, 
containing the remaining 24 free parameters. For fitting the neutrino data, however, it is 
more useful to first rewrite the neutrino Yukawa couplings using a generalization of the 



Casas-Ibarra parametrization 



64] 



Consider first the effective mass matrix of the light neutrinos for the inverse seesaw. It 
is given by 



m 



eff 



We can rewrite m D as 



m D = M^V; T (V^)" 1 J RDiag( v ^v)^. 



(16) 



(17) 



Here U v is the mixing matrix determined by the oscillation experiments, in the basis where 
the charged lepton mass matrix is diagonal, m Vi are the three light neutrino masses, R is an 
arbitrary real orthogonal 3x3 matrix and fis are the eigenvalues of the matrix /i^ with 
the matrix which diagonalizes fis- 

Eqs ( 1T61) and ( jTTI) allow to fit neutrino data in a straightforward way, if the tree-level 
contribution dominates, see below. Since one can always choose a basis where Mr is diagonal, 
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the flavour violation necessary to fit oscillation data resides in mj) and in fis- Particularly 
simple solutions are found, assuming either /ij or mo are diagonal too. For diagonal [As, for 
example, one finds 

— (lg) 



m D = Diag( 



■M Ri )U v . 



Oscillation experiments have determined the mass squared differences and mixing angles 
of the active neutrinos with high precision, see for example js]. Recently also the last of the 
mixing angles in the left-handed neutrino sector has been measured in two reactor neutrino 



experiments, DAYA-BAY 



66] and RENO 67J. With all these data, the situation can be 



summarized as follows: The atmospheric neutrino mass squared difference and angles are 
A(mi tm ) = (2.31 - 2.74) x 1(T 3 eV 2 (normal hierarchy) and sin 2 # A tm = 0.36 - 0.68, the 
solar parameters are A(m|) = (7.12-8.20) x 10~ 5 eV 2 and sin 2 0© = 0.27-0.37 and finally 
sin 2 6*i3 = 0.017 — 0.033, all at 3 a c.l. js]. Apart from the data on the reactor angle, neutrino 
angles are still well-fitted by the tribimaximal mixing ansatz {ssj , which has sin 2 #Atm = 1/2 
and sin 2 (9© = 1/3. 

The large atmospheric and solar angles require large off-diagonals in at least one of the 
two matrices Y v or For the case of strict normal hierarchy (m vi = 0) and diagonal fi$, 
oscillation data can be well fitted to leading order in the small parameter sin # 13 by 

/ n n n \ 



Y = \Y 





a 

V^sinfl 





a(i-^; 



13 





-a(l + 
1 



sin 6>i3 ' 
V2 - 



(19) 



with 



a= (Aml/Ami tm )* -0.4 



(20) 



where \Y V \ can be easily calculated from fi$ and Mr. 

The above discussion is valid at tree-level. In the inverse seesaw neutrino masses also 



receive important corrections at the 1-loop level, once B fls becomes sizeable 69) . An example 
is shown in fig. (JT]). Here, we have chosen as an example mo = 100, Mi/ 2 = 1000, Aq = 
(all in GeV) tan/3 = 10, sgn(/x) > and Mr = 250 GeV. For this plot we assume fig and 
Bm r = 3 x 10 4 GeV 2 to be diagonal and degenerate. Y v is then fitted by eq.f ll7l) . A smaller 
value of fis implies then a larger value for the entries in Y v . 

In the left of fig. (pQ) we show m V2 and m V3 as function of -B Mg , while the plot on the 
right shows the same neutrino masses as a function of m 2 ^, the mass squared difference 
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US [GeV]: 4.2 x 10" 10 5.6 x lCT' 1 5.6 x 1CT 7 



< 




/is [GeV]: 4.2 x 10" ln 5.6 x 10" 9 5.6 x 10" 7 



> 




10" 3 10 



IQ-' VT' W 10" 

m/, [GeV] 



FIG. 1: Neutrino masses versus B^ s (left) and versus (right) for one particular but arbitrary 
parameter point (see text), for three different values of fj,s- 

between the CP-even and CP-odd sneutrinos. This splitting is proportional to -B MS , while 
loop neutrino masses are proportional to Y^B^ S . To restrict the neutrino mass to be smaller 
than the atmospheric scale than results in an upper limit on Y^B^g. For \\is\ ~ 5 x 10~ 7 
GeV, corresponding to the largest entries in Y v to be of order O(10~ 2 ), the splitting can be 
as large as C(10 _1 ) GeV. Note, however, that with typical mSugra-like boundary conditions 
one expects naively that B^ s ~ Hs m o ~ 10~ 4 — 10~ 7 GeV 2 . In this case splitting between 
the sneutrinos becomes negligible. 

2. Linear seesaw 



For the linear seesaw the neutrino mass matrix is given by 

/ m T D Ml\ 
M v = m D M R 
\M L Ml / 

with the effective neutrino mass matrix for the light neutrinos given as 



m v = m T D M T R 1 Mi + MjM^ l m D . 



For the linear seesaw one finds for the CI parametrization 65] 



(21) 



(22) 



(23) 
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where A has the following general form: 



( 1 a 

i 



V 



—a 2 c 



(24) 



with a, b, c real, but arbitrary numbers. The parameter counting for the linear seesaw is 
given in table ( 1IIII) . We have in total 36 real parameters, 24 moduli (12 masses and 12 
mixing angles) plus 12 phases. Fits to neutrino data can be easily done using eqs (122]) and 



For example, for strict normal hierarchy, degenerate Mr and diagonal and degenerate 
Ysl one finds to leading order in #13 

/ 



\Yu\{ 







sin 6>i3 sin6>i3 ^ 

sin 0i3 1 1 
2 2 



V2 
sin 6>i3 
V V2 

( 



(25) 



/ 



m 

H x 

3 



13 X 



u 2 2 

\ 1 



)>• 



-i-i i \ 

-1-11 + v^sin^ 

v 1 1 -V 

where again, the prefactor \Y U \ can be calculated from \Ysl\ and Mr. Note that the flavour 
structure of eq. ( 1251) differs significantly from eq. ( TT9l) for the same choice of angles, see the 
discussion about lepton flavour violation in the next subsection. 





Parameters 


Moduli 


Phases 


Ny 


Y e ,Y u ,Y SL ,M R Diractype 


4 x n 2 


4 x n 2 


G 


U(n) L U(n) u c [/(n) e c (g) ^7(n) 5 


4 x n(n - 1} 


4 x n(n 2 +1) 


G' 


no LF conservation 






Nphys 




24 
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TABLE III: Parameter counting for the linear seesaw model for three generations. 



In case of the linear seesaw, loop contributions to the neutrino masses from the splitting in 
the sneutrino sector is always negligible for neutrino masses in the sub-eV range, assuming 
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FIG. 2: Br(^ — > ej) calculated in the inverse (left) and linear (right) seesaw. Here, all flavour has 
)een put into the Yukawa Y u , while neutrino angles have been fitted to their best fit point values 
8j. A random scan over mo and M R in the interval [100,1000] GeV and B Mr = [10 3 , 10 6 ] GeV 2 
has been performed for these points. Note that the axis are the same for inverse and linear seesaw 
for an easier comparison. Linear seesaw leads to smaller LFV than inverse seesaw for equal choice 
of neutrino angles. 



the trilinears to be proportional to T x oc Y X A . This can be understood as follows: The 
difference in the eigenvalues of the CP-even and CP-odd sector is entirely due to the different 
signs in the off-diagonals in eq. (Q. As can be easily shown, the maximum difference in 
the eigenvalues is reached for Y S l — Y v . However, eq.( )22|) shows that the product Y SL Y U is 
required to be small, due to the observed smallness of neutrino masses. Thus, the splitting 
in the sneutrino sector in case of linear seesaw is maximally of the order of m u msusY, i-e. 
0(lCr 9 ) GeV 2 . 



B. Lepton flavour violation 



In any supersymmetric model, limits on lepton flavour violating decays such as fi — > 
provide an important constraint on the parameter space [tt1]- In models with a low scale 
seesaw especially important constraints come from U — > 3lj 



in nuclei 



72|. 



7l| and from \i — e conversion 



The fit to neutrino data requires non-trivial flavour violating entries in at least one of the 
Yukawa- or mass matrices: Y v or Y$l for linear and Y v or \xs for inverse seesaw. If we assume 
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that the LFV resides in Y v , limits on the Yukawa result as shown in fig. (J2J). In this figure 
we have chosen /j,s (left) or Y$l (right) diagonal and neutrino angles have been fitted to their 
best fit point values J§] using Y u . A random scan over m and Mr in the interval [100,1000] 
GeV and Bm r = [10 3 , 10 6 ] GeV 2 has been performed for these points, fixing tan /3 = 10 and 
M1/2 — 2.5 TeV. Upper limits of the order of (few) 10~ 2 (10 _1 ) result for inverse (linear) 
seesaw, despite the heavy SUSY spectrum (due to the large value of M1/2). Much stronger 
limits result for lighter spectra. Note that U — >■ 31 j |7l| and /jl — e conversion in nuclei [72 1 
can lead to even stronger limits. We will not repeat this exercise here. 

Note also, as discussed in the next section, that the constraints from relic density of 
sneutrinos lead to an approximate lower bound on the absolute size of the Yukawa coupling 
\Y± 



IV. SNEUTRINO DARK MATTER 



In this section we discuss the relic abundance (RA) and the direct detection cross section 
(DD) of sneutrinos in the different models. We will first discuss the simpler case of the 
inverse/linear seesaw and then turn to the mBLR model. 

In order to reduce the number of free parameters in our numerical scans, we calcu- 
late all spectra with CMSSM-like boundary conditions, i.e. at the GUT scale we choose 
(m , M1/2, A)> tan/3, sgn(/i)), from which all soft parameters at the electro-weak scale are 
calculated using full 2-loop RGEs. Unless noted otherwise, we always assume that the tri- 
linear soft parameters are related to the superpotential parameters in a "mSugra"-like way: 
T a oc Y a A at m GUT . 

In addition to the MSSM parameters, we have the neutrino Yukawa couplings Y u and 
several model specific parameters. These are Mr and B Mr and, in case of the inverse 
(linear) seesaw li s and B^ s (Y SL ). While, in principle, all of these are matrices we use eq. 
f fT7|) and ( 123]) ) to fit neutrino data and usually assume all matrices are diagonal except one. 

For the mBLR model we have the free parameters Y s , vr, tan/3fj, fi R and rriA R - Recall, 
Mr = Y s vr/\/2 and vtia r is the CP-odd scalar Higgs mass in the \r sector. Due to the 
constraints from LFV discussed above, we usually put all LFV into either fi$ (inverse seesaw) 
or Ysl (linear seesaw). This way we only have to check for the constraints from Z° and Higgs 
physics and lower limits on squarks and gluinos discussed in section (JXTTJ) . 
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A. Inverse/Linear seesaw 



Sneutrinos can be the LSP, practically independent from the actual choice of the CMSSM 
parameters. This can be easily understood from eqs ([7]) and (jU]) and is demonstrated by two 
simple examples in fig. flS]). 

In fig. 03]) we show two examples of tree-level sneutrino masses calculated as function of 
Bm r for two particular but arbitrary choices of parameters: mo = 100, Mu 2 = 1000, A = 
and // = 800 all in GeV and \Y V \ = 0.1 and tan/3 = 10. In addition, M R = 200 GeV (left) 
and Mr = 500 GeV (right). This calculation was made in a one generation toy model and 
serves only for illustration. The general behaviour is easily understood. First, recall that 
within CMSSM roughly m x o ~ ~ 0AMi/ 2 - Entries on the diagonals of the sneutrino 
mass matrix are of the order m\ L ~ m\ + 0.5Mjy 2 , vr? v c v c — m\ + M R and m| 5 ~ m\ + M R . 
If a/ m\ + M\ < 0.4Mi/ 2 (one of the pair of) right sneutrinos is the LSP, see left plot. On 
the other hand, for larger values of mo and or Mr, right sneutrinos still can be the LSP 
if Bm r }t \Z m o + M R , since in this case a large off-diagonal in the sneutrino mass matrix 
leads to a large splitting between the two lightest eigenstates, with the lighter one becoming 
very light, see right plot. Since Bm r is naively expected to be of order m 2 SUSY) sizeable 
splitting between the right sneutrinos is expected and in a random scan over parameters 
such sneutrinos emerge as LSP quite often. Note, that a light eigenvalue in the sneutrino 
sector can also be made by a large off-diagonal in the sneutrino mass matrix in the LR and 
LS entries of the mass matrix. 

In the early universe sneutrinos can annihilate into SM particles through various types of 
interactions. The most important Feynman diagrams are shown in figs (jl]) and ([5]). Fig. (j4j) 
shows the quartic interaction between two sneutrinos and two Higgses and s-channel Higgs 
exchange. The former is very efficient for m„ LSP > m h o, while the latter is important near 
mo LSP — m h o/2. Fig. ([5]) shows the quartic interaction with W- and Z-bosons and t-channel 
neutralino exchange. The importance of the latter depends on the SUSY spectrum. 

The relative importance of different diagrams is strongly dependent on the kinematical 
regime. A typical example of final state branching ratios versus the lightest sneutrino mass 
is shown in fig ([6]). In this scan we have fixed m = 120, M\m = 600, A = all in GeV and 
\Y„\ = 0.4 and tan/3 = 10. In addition /i 5 = [10 -11 , 10~ 9 ] GeV. 

From left to right we see that the most important channels are vlsp^lsp — ^ T ^ (ma- 
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FIG. 3: Two examples of tree-level sneutrino masses calculated as function of Bm r for two partic- 
ular but arbitrary choices of parameters: mo = 100, M]y 2 = 1000, Aq = and [i = 800 all in GeV 
and \Y V \ = 0.1 and tan/3 = 10. In addition M R = 200 GeV (left) and M R = 500 GeV (right). For 
comparison also the lightest neutralino mass is shown. 
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FIG. 4: Examples of Feynman diagrams contributing to the vlsp^lsp annihilation: To the left 
quartic interaction; to the right s-channel Higgs exchange. 

genta with triangles), v LS pv LS p — > bb (brown), v lsp v L sp — > vv R (orange), v LS pVlsp — > 
W + W~ (red ), v LSP v LSP — > Z°Z° (purple), v LSP v LSP — > HH (blue), v LSP v LSP — > tt 
(green); finally, to the right of the figure, the contributions coming from the coannihilations 
are shown: ee — > tt (magenta with triangles), and ee — > 77 (in yellow). 

For low sneutrino masses the determination of the relic abundance is dominated by Higgs 
exchange, see fig. (HI) right. Since the Higgs couplings are proportional to SM fermion 
masses, bb is most important in the low mass regime, followed by tt. For sneutrino masses 
above approximately rn„ LSP ~ 45 GeV the final state vv becomes dominant in this example. 
This is because with these parameter choices the lightest of the "singlet" neutrinos has a 



17 




FIG. 5: Examples of Feynman diagrams contributing to the vlsp^lsp annihilation: To the left 
quartic interaction with gauge bosons; to the right t-channel neutralino exchange. 
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FIG. 6: Examples of final state branching ratios for the annihilation cross section of sneutrinos to 
SM final states versus the lightest sneutrino mass (in GeV). For the parameter choices of this scan, 
see text. Calculation uses the inverse seesaw model. Different kinematical regimes are visible, see 
discussion. 

mass of about 45 GeV and the Higgs couples always to vl^r-, he. one light and one heavy 
neutrino. 

Single Z° exchange is less important than Higgs exchange, since scalar-scalar-vector cou- 
plings are momentum suppressed. For m„ LSP > 80 GeV, however, two gauge boson final 
states become dominant, the channel W + W~ being more important than Z°Z°. For masses 
above vny hSP > 120 GeV also two Higgs final states are sizable. All these final states are due 
to quartic interactions, see fig. (BJ left and fig. ([5]). Due to the large top Yukawa coupling, 
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FIG. 7: General scan for the inverse seesaw model. The plot shows O/i 2 versus the mass of the 
lightest sneutrino (in GeV) for points in which the sneutrino is the LSP. 

the two top final state, once kinematically possible, becomes very important. And, finally, 
for rriv LSP approaching the NLSP mass, in this example the lightest scalar tau, coannihilation 
into taus becomes dominant. 

Next, we have performed a general scan over the parameter space of the model choosing 
randomly (m , M 1/2 , A , tan/3, sgn(/i)) in the interval m = [100,3000], M 1/2 = [200,3000], 
A = 0, tan/3 = 10 and sgn(/i) > and \Y V \ = 0.3, M R = [0,1000]. B Mr is calculated 
accordingly to enhance the percentage of sneutrino LSP points. We post-select data points 
with sneutrino LSPs and cut on all points not fulfilling the lower bounds on squark and 



gluinos masses from the LHC 61J . Results are shown in fig. ((Tj) for the case of the inverse 
seesaw. Shown is the calculated RA (Qh 2 ) versus the mass of the lightest sneutrino for points 

n 

in which the sneutrino is the LSP. The band, which is the allowed range from WMAP [731 ]. 
shows that one can easily get points with the correct relic abundance over a wide range of 
parameters. The figure is for the inverse seesaw, linear seesaw is qualitatively very similar. 

The plot shows several distinct features. First, for masses of sneutrinos around my LSP ~ 
60 GeV a strong reduction of the RA occurs, due to the s-channel Higgs exchange. As can 
be seen, this diagram is very effective in reducing the RA whenever m^ LSP is within a few 
GeV of the mass of the Higgs, but less important elsewhere. In the region above my LSP = 80 
GeV, quartic interactions with the gauge bosons are effective and above mc LSP = 175 GeV 
two-top final states become dominant. For very large Tnp LSP one sees an overall trend that 
the RA rises with rising sneutrino mass, apart from a few scattered points. Low RA, i.e. 

19 




10" 6 10" 5 10" 4 10" 3 10" 2 10"' 10° 10" 6 10" 5 10" 4 10~ 3 10" 2 10"' 10' 



\Y \ \Y \ 

FIG. 8: Scan for the inverse (left) and linear (right) seesaw model. The plot shows Qh 2 versus 
\Y U \ in a scan over the remaining parameters, see text. The color coding of the points shows the 
mass difference between the lightest sneutrino mass and the NLSP (next-to-LSP) mass, in this scan 
practically always the lightest of the charged sleptons. Red: uinlsp — i^lsp < 100 GeV, violet: 
100 < rriNLSp — itilsp < 500 GeV, blue: tunlsp — m<LSP > 500 GeV. For a discussion see text. 

Qh 2 ~ 0.1, in this high mass regime can practically only be made via co- annihilation or 
s-channel heavy Higgs exchange. Note that the fact that there are only a few points with 
m i>LSP below 50 GeV is just an artifact of the scanning procedure. However, the general 
trend is that for very light sneutrinos the calculated RA is larger than Qh 2 ~ 0.1. We will 
come back to a more detailed discussion of light sneutrinos in the next section. 

In fig. ([Tj) we have fixed the neutrino Yukawa couplings to a constant value. However, 
sneutrinos which are purely singlets do not couple to gauge bosons and thus their relic 
abundance is usually too large. For mixed sneutrinos the RA depends strongly on the 
choice of \Y U \. An example is shown in fig. (jSJ. The figure shows on the left (right) results 
for the inverse (linear) seesaw. In both cases we have fitted neutrino data, using eqs (JT7J) 
and (1231) . and scanned over the parameters: mo and Mr in the interval [100,1000] GeV and 
B Mr = [10 3 , 10 6 ] GeV 2 . Here, M 1/2 was fixed to M 1/2 = 2.5 TeV and tan/3 = 10 and A = 0. 
The choice of such a large M\i2 guarantees that all points have a lightest Higgs mass in the 
vicinity of 125 GeV. It also makes all SUSY particles, except the sneutrino, relatively heavy. 

The points in fig. (jHJ) are color coded by the mass difference between the lightest sneutrino 
mass and the NLSP (next-to-LSP) mass, in this scan practically always the lightest of the 
charged sleptons. Red: m NLSP — m LS p < 100 GeV, violet: 100 < m^LSP ~ m LSP < 500 
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FIG. 9: Diagrams contributing to the direct detection cross section: elastic scattering of vlsp ° ver 
quarks. 

GeV, blue: ttinlsp ~ m LSP > 500 GeV. For large \Y U \ the RA goes down as flh 2 oc |1^|~ 4 , 
for small values of Y v the points show practically no dependence on \Y U \. This is because 
the determination of the RA is then dominated by coannihilation processes with the lightest 
stau. These can be very efficient, if Am 2 = m 2 — vn^ hSP 2 ~ few GeV, less so for larger mass 
differences. Thus, to reduce the relic density of the sneutrino to acceptably small values, one 
needs either a special kinematic configuration, such as co-annihilation or s-channel resonance, 
or \Y V \ has to be larger than roughly \Y V \ > 0.1. 



1. Direct Detection 



Direct detection of the sneutrinos consists in detecting the recoil energy coming from the 
elastic scattering of sneutrinos with nuclei inside a detector. The interaction, which occurs 
in the non relativistic limit, since the velocity of dark matter particles in the Galactic halo is 
small, comes from basically two diagrams contributing at tree level: the t-channel exchange 
of a neutral Higgs or of the Z boson. Which of the two diagrams is the more important one 
depends on the actual value of \Y V \. 



The Z-boson exchange cross section is 30] : 



f ^ir+tp ^ + 2(2sii ^ - i)z ^ (2e) 

where is the nucleus mass, A^- and are the mass number and proton number of the 

21 



nucleus and Gf is the Fermi constant. 

The Higgs-bosons exchange scattering cross section is |30J 



a 



Higgs 
vlsp A/" 
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VLSP 



[fpZtf + fn {AjKf - Ztf)Y 



(27) 



where Af denotes the nucleus, and the quantities and Zjj are the mass number and 
proton number of the nucleus, f p and /„ are hadronic matrix elements which parametrize 
the quark composition of the proton and the neutron, and which represent the effective 
coupling of the v LS p to the nucleus, but are subject to considerable uncertanties 3(], 7o1 |. 
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FIG. 10: Direct detection cross section (in [cm 2 ]) for sneutrino LSPs (masses in GeV), for the 
inverse seesaw model. The points are those from fig. ([7]) compatible with the upper bound on the 
relic abundance. Also the current limits from XENON-100 [7_J] (red line), CDMS [76[ (green line), 
DAMA (with and without channeling, orange regions) [^j], and Cogent [78] (purple region) are 
shown for comparison. 

In figure ( flOl we depict the direct detection cross section versus the LSP sneutrino mass 
(blue points). The points are the same as shown in fig. ([7]), but after a cut on the relic 
abundance. In the same plot, the current limits from XENON-100 |74| (red line), CDMS 
76 1 (green line), DAMA (with and without channeling, orange regions) [tT ] , and Cogent 
78| (purple region) are shown. The major bound nowadays comes from the XENON-100 



experiment 



74( , whose best sensitivity is around 10 44 cm 2 for a dark matter candidate of 50 



u LSP > ioo 



GeV. The sneutrinos show a SI cross section asj < 10 _42 cm 2 , and for masses m 
GeV they are compatible with current limits by XENON-100. However, XENON-IT, whose 
sensibility should improve up to 10~ 46 cm 2 , will test those cross sections. 
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We have not been able to find low sneutrino masses of the order of (9(5 — 10) GeV, 



which have the correct relic density and fulfi 



direct detection experiment XENON-100 74|. However, this calculation has been done 



1 at the same time the constraints from the 



with B^g oc rriQfis and lepton number violation in the sneutrino mass matrix leads to the 
mass splitting between the real and the imaginary part of the lightest sneutrino, and the 
scattering via Z boson exchange occurs inelastically, through a transition from the real to 
the imaginary or viceversa. Points shown in fig. ([TO]) have all very small splitting in the 
sneutrino sector, but if the mass splitting is greater than some keV, scattering is strongly 
suppressed at direct detection experiments. Indeed, the maximum kinetic energy that the 
sneutrino LSP can transfer to the detector depends on the velocity it moves relative to the 
nucleus v (~ 10~ 3 in the galactic halo), the nucleus mass M and the angle 6 of scattering: 

A 2 2 

E = ^(l-cos(9)) (28) 

where A = mvLSpM ,, , which would give, in the case of a Xenon detector for instance, and 
m i>LSP = 100 GeV, E = 25 keV (if cos(#) = 0). For heavier sneutrinos with a mass of the 
order of TeV, for a splitting larger than some hundred keV the direct detection cross section 
goes to zero. Such "large" splitting is currently not excluded in the inverse seesaw, compare 
to fig. ([TJ). Thus, in principle inverse seesaw can evade all constraints from direct detection, 
while linear seesaw can not, see the discussion in section ( IIII A2j) . 



B. mBLR model 



In this subsection we discuss the DM phenomenology of the supersymmetric U(1)r x 
U(1)b-l extension of the standard model. The main difference to the simpler models dis- 
cussed previously are the presence of the extra gauge boson Z' and the possibility to have 
an additional light, mostly singlet Higgs boson, which lead to some important changes in 
the phenomenology. 

First, recall that the U(1)r x U(1)b-l gauge symmetry of this model is spontaneously 
broken to the hypercharge group U(l)y by the vevs v XR and v% R of the scalar components of 
the xr and xr superfields whereas the SU(2)l <8) U(l)y — > U(1)q breaking is driven by the 
vevs Vd and v u of the neutral scalar components of the SU (2) l Higgs doublets Hd and H u up 
to gauge kinetic mixing effects. The tadpole equations for the different vevs can be solved 
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FIG. 11: Scan for the CmBLR version of the extended gauge model. Parameters are varied as 
follows: m = [0,6000] GeV, M 1/2 = [3000,8000] GeV, tan/3 fl = [1.0,1.3]. The other parameters 
are set to the values tan/3=10, Aq = —4500 GeV, Ys = diag(0.3); vr has been chosen different in 
the two plots, vr = 6 TeV and vr = 10 TeV, respectively. Masses of the vlsp are in GeV. 



analytically for either (i) (fj,, B^, /j,r, B m ) or (ii) (//, B^, m 
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We address the minimal version option (i) as CmBLR (constrained mBLR), since it allows 
to define boundary conditions for all scalar soft masses at mcuT-, reducing the number of free 
parameters by four, although leading to some constraints on the parameter space, such as a 
lower bound on > 1) [TJ]. The second option (ii) is instead more flexible, and 

we have made use of it in some of our scans, too. We will refer to this option as x/? m BLR 
version (non- universal xr masses mBLR). We have not used the last option, which we only 
mentioned for the sake of completeness. 

The result of Qh 2 for two general scans is shown in fig. ffTTT) . Parameters have been 
scanned as described in the figure caption. Note that there are two fixed but different 
choices of vr in the left and right plots, leading to two different values of the Z' mass. In 
both plots in fig. flllD the main feature clearly visible is the Z' pole. Indeed, the annihilation 
of the vlsp LSPs into SM particles via the Z' becomes efficient when the mass of the Vlsp is 
close to half the mass of the Z' . The mass of the Z' can be calculated analytically [13] and 
mainly depends on the value of vr, see eq. (|T3|) . The ATLAS searches for a Z' set a lower 
limit on its mass which is 1.8 TeV, and this translates into a lower limit on v r > 5 TeV, see 
the plot on the left. The plot on the right shows that choosing a higher value of vr we can 
get very heavy vlsp DM with the correct RA, up to masses of several TeV. 
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FIG. 12: Final state branching ratios for the annihilation cross section of sneutrinos to SM final 
states versus the lightest sneutrino mass (in GeV). For the parameter choices of these scans see fig. 
pip . For a discussion of the different kinematical regimes which are visible, see text. 

The main annihilation channels for sneutrino DM for the points of fig. (TTTj) are shown 
in fig. f lT2|) . Far from the Z'-pole resonance these are vlsp^lsp — * T ^ (magenta ), 
vlspVlsp — ► bb (brown), v L sp^lsp — > vv R (orange), v LSP v LSP — > W + W~ (red), 
vlspVlsp — > Z°Z°/Z R Z R (purple ), vlsp^lsp — > HH (blue), v lsp vlsp — > it (green). 
The quartic coupling with two Higgses ( h°, h° BLR and A , depending on if they are kine- 
matically allowed, depending on the vlsp mass) is one of the most efficients, as before. For 
lower masses the annihilation via the MSSM Higgs is the most efficient, as can be noticed 
by the small relic density for lower masses, expecially in the first plot, where on the left 
end side we are approaching the region where the quartic Higgs coupling is important (for 



m 



vlsp 



120 GeV). 



Recall that in this model the Higgs sector is more complicated due to the extended gauge 
structure. The U(1)b-l x U(1)r breaking results in one additional light Higgs, h° BLR |l2 |. 
The mixing between the MSSM Higgs h° and the h° BLR enhances the mass of the mostly 
MSSM Higgs, leading to a MSSM-like Higgs in accord with the most recent ATLAS and 
CMS preferred regions, without much constraints on the SUSY spectrum. However, this 
enhancement of the MSSM Higgs mass occurs usually in the model if the h° BLR has a mass 
of the order of the MSSM-like state or less, i.e. the presence of a light singlet Higgs is 
preferred unless the SUSY spectrum is rather heavy (in which case the CMSSM limit is 
reached) . 

We have also checked for constraints coming from direct detection in the limit of negligible 
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FIG. 13: Direct detection cross section (in cm 2 ) for sneutrino LSPs in the BLR model. Masses 
are in GeV. Black points refer to the scan described in fig. (jlip left with vr = 6 TeV. Blue points 
stand for the scan of fig. (|lip right with vr ; = 10 TeV. All points shown fulfill the RA constraints. 
Higher vr leads in general to lower DD cross section. 

sneutrino splitting. Examples for direct detection cross section are shown in fig. ( TT3l) . As 
before, see fig. fllOp . different experimental constraint are also shown. All points shown fulfill 
the constraints from relic abundance. We have calculated two scans, one with vr = 6 TeV 
(black) and one with Vr = 10 TeV (blue). As can be seen, practically all of the points with 
Vr — 6 TeV in this scan are excluded by the limit from XENON-100, while most of the 
Vr = 10 TeV are allowed. Thus XENON-100 puts currently a lower bound on Vr (and thus 
the Z' mass) of the order of vr ~ 10 TeV for sneutrino LSPs as DM. 

The origin of this surprisingly strong constraints lies in the Z° — Z' mixing. The mixing 
angle between these two states is roughly of the order 9 z oz> ~ (9lv 2 ) / \9r v r)- Thus the 
right sneutrinos, which couple mostly to the Z\ couple via this mixing also to the Z° . The 
Z° has an experimentally fixed mass. Thus, the only possibility to suppress the DD cross 
section 4 is to increase v r. 

Finally, fig. f[T4"j) shows a dedicated scan for low mass sneutrinos in the mBLR model. The 
different curves are slight variations of the parameters near the study point BLRSP1. The 
original parameters of BLRSP1 were: m = 470, M l/2 = 700, tan/3 = 20, A = }Vr = 4700, 

4 Apart from a large sneutrino splitting. 
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FIG. 14: Scan into the low sneutrino mass region using the mBLR model. For the parameter 
choices see text. The figure demonstrates that the mBLR mode can give the correct RA for low 
mass sneutrinos in those parts of the parameter space where a light, singlet Higgs is present. 

tan/?R = 1.05, lir = —1650 and Ma r = 4800 GeV. To obtain very light sneutrinos, m has 
been lowered to mo = 440 GeV, while the different curves are for Mi/ 2 = 650, 660, 675 and 
700 GeV and the scanning runs vn^ R from 3000 — 4000 GeV. The resulting scan produces 
sneutrinos with masses in the interval [5, 100] GeV, while the lightest Higgs mass, in this case 
a mostly singlet Higgs, has a mass eigenvalues of — 1 — 50 GeV. The figure shows a pole 
around m„ LSP ~ 62 GeV, due to a mass for the MSSM-like Higgs of around 124 — 125 GeV 
in all cases. There appear additional dips in the RA for smaller sneutrino masses, whenever 
m^ LSP ~ m^/2. This demonstrates that in the extended gauge model it is possible to have 
the correct RA even for very low sneutrino masses. 

However, note that, while the model can in principle give DD cross section large enough 
to explain the DAMA [tT ] , and Cogent Es] hints, such points will always be inconsistent 



with the constraints from XENON-100 [74]], also for the case of inelastic dark matter 79] 



V. CONCLUSIONS 



We have studied low scale seesaw models with a sneutrino LSP. We considered two pos- 
sibilities: Models with the MSSM gauge group and either a linear or inverse seesaw and a 
model with the gauge group SU(3) C x SU(2)l x U(1)b-l x U(1)r and an inverse seesaw. 
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Sneutrinos can be the DM in both cases, fulfilling all known experimental bounds. 

However, while inverse and linear seesaw lead to different results for LFV, in general, 
they give similar DM results. There are some differences in detail, though: In the inverse 
seesaw it is possible to avoid all direct detection constraints using a large enough splitting 
in the sneutrino sector, which leads to "inelastic" dark matter. This is not possible in the 
linear seesaw, due to constraints from neutrino physics. 

In the extended gauge model there is more freedom than in the simpler MSSM-group 
based models. Especially very light (0(1) GeV) or very heavy (O(several) TeV) sneutrinos 
can give the correct relic density, due to the existence of a mostly singlet Higgs in the former 
case and due to the Z' in the latter. Very light sneutrinos could explain the hints from 



DAMA 



77j or Cogent (78|, but are inconsistent then with XENON-100 [7J, |79| 



Finally, it is interesting to note that in the limit of small sneutrino mass splitting the 
DD limit from XENON-100 Q leads to a lower limit on v R of the order of O(10) TeV for 
sneutrino LSPs as the dominant component of the galactic dark matter. 
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